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Abstract. We establish long-time stability of multi-dimensional viscous shocks of a gen- 
eral class of symmetric hyperbolic-parabolic systems with variable multiplicities, notably 
including the equations of compressible magnetohydrodynamics (MHD) in dimensions 
d > 2. This extends the existing result established by K. Zumbrun for systems with 
characteristics of constant multiplicity to the ones with variable multiplicity, yielding the 
first such a stability result for (fast) MHD shocks. At the same time, we are able to drop 
a technical assumption on structure of the so-called glancing set that was necessarily used 
in previous analyses. The key idea to the improvements is to introduce a new simple 
argument for obtaining a L 1 — > L p resolvent bound in low-frequency regimes by employ- 
ing the recent construction of degenerate Kreiss' symmetrizers by O. Gues, G. Metivier, 
M. Williams, and K. Zumbrun. Thus, at the low-frequency resolvent bound level, our 
analysis gives an alternative to the earlier pointwise Green's function approach of K. Zum- 
brun. High-frequency solution operator bounds have been previously established entirely 
by nonlinear energy estimates. 
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1. Introduction 
We consider a general system of viscous conservation laws (d > 2) 

(1.1) U t + J2F j (U) Xj =^2( B3k (U)U Xk ) X] , xeR d , t>0, 

3 jk 

U,F* G R n , B jk € R nxn , n > 2, with initial data U(x,0) = U (x), and a planar viscous 
shock, connecting the endstates U±: 

(1.2) U = U{x 1 ), lim U{x 1 ) = U±. 

We study the long-time linearized and nonlinear stability of the viscous shock U un- 
der multi-dimensional perturbations of initial data. The problem has been carefully and 
successfully investigated by K. Zumbrun and his collaborators in \Z2\ \Z3\ \Z4\ IGMWZl] . 
There, due to technical arguments of the analysis, the authors put assumptions on the mul- 
tiplicity of hyperbolic characteristic roots and structure of the so-called glancing set (see 
(H4)-(H5) below). The latter assumption (which is automatically satisfied in dimensions 
d = 1,2 and in any dimension for rotationally invariant problems) assures the glancing set 
to be confined to a finite union of smooth curves on which the branching eigenvalue has 
constant multiplicity. This is precisely to reduce the complexity of multi-variable matrix 
perturbation problem when dealing with glancing blocks to a simplified form of a two- 
variable perturbation problem. Whereas, the constant multiplicity assumption excludes an 
important physical application, namely, the equations of magnetohydro dynamics (MHD) 
in dimensions d > 2. In the current paper, we are able to relax the assumption of constant 
multiplicities to variable multiplicities, allowing (fast) MHD shocks to be treated and thus 
yielding for the first time the long-time multi-dimensional stability for these shocks. In 
addition, we are also able to drop the assumption on structure of the glancing set at a price 
of having t 1//4 slower in decay rates in dimensions d > 3. 

Our main improvements rely on recent remarkable and technical works of O. Gues, G. 
Metivier, M. Williams, and K. Zumbrun [GMWZ5, GMWZ6J where the authors have ob- 
tained the 1? stability estimates and small viscosity stability for the symmetric systems 
with variable multiplicities via their construction of Kreiss' symmetrizers. The idea is to 
employ these available estimates to establish the long-time stability, or more precisely, to 
derive a resolvent bound in low-frequency regimes. This will be the main contribution of our 
present paper. High-frequency estimates are already established by K. Zumbrun via elegant 
nonlinear energy estimates for a very general class of symmetrizable systems, including our 
class under consideration. 



STABILITY OF MULTI-DIMENSIONAL VISCOUS SHOCKS 



3 



We would like to mention that the idea of using L 2 stability estimates via the construc- 
tion of degenerate Kreiss' symmetrizers to attack the long-time stability problem has been 
investigated in [GMWZlj. There the authors obtain the result under (H4)-(H5) assump- 
tions (and treat the strictly parabolic systems). In our analysis, we avoid these technical 
assumptions, by introducing a rather simpler argument for L 1 — > LP resolvent bounds in 
low-frequency regimes, which turns out to be the key to the improvements. The analysis 
works precisely for the case of dimensions d > 3. In dimension d = 2 (the condition (H5) 
is now always satisfied), the analysis of [GMWZ1] indeed works even for the MHD shocks 
as we are considering here by combining their later work in [GMWZ6] (though it was not 
stated there). In Section [H we represent a slightly modified version of [GMWZ1] treat- 
ing this two-dimensional case, or more generally, cases with (H5) in a more direct way. 
Once these low-frequency resolvent bounds are obtained, the stability analysis follows in a 
standard fashion \L2\ IZ31 [Z4], See Section [L6l for further discussions. 

1.1. Equations and assumptions. We consider the general hyperbolic-parabolic system 
of conservation laws (jl.ip in conserved variable U, with 

v=(#i)> Bik= (bf &?< 

u 1 E R n - r , u 11 € R r , and 

(1.3) ^^4%£fc>0|£l 2 >o, V£eM n \{o}. 

jk 

Following [Z3|. IZ4] . we assume that equations (11. ip can be written, alternatively, after a 
triangular change of coordinates 

(1.4) W := W(U) = ( ^f'flj,) , 

\w (it , u )) 

in the quasilinear, partially symmetric hyperbolic-parabolic form 

(1.5) A°W t + AiW Xj = J2(& k W Xk ) Xj + G, 

where 

20 _ \ ~ x _ (A\ x A\ 2 \ ~ jk _ (0 

\ A\ 2 ) ' " \A\ X AlJ ' " V ^ 

and, defining W± := W(U±), 

(Al) A^(W ± ),A°,A\ l are symmetric, A > 6 > 0, 

(A2) for each ( £ 1^ {0}, no eigenvector of Ylj t,jA° ' {A°)~ l (W±) lies in the kernel of 
E jk ^k& k (A»)-i(W±), 

(A3) SfcrEfr'*^* > 8\H\ 2 , and G = (?\ with g(W x ,W x ) = <D(\W X \ 2 ). 

Along with the above structural assumptions, we make the following technical hypotheses: 
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(HO) Fj,Bi k ,A°,Ai,Bi k ,W(-),g(-,-) G C s+1 , for s > [(d - l)/2] + 2 in our analysis of 
linearized stability, and s > := [(d — l)/2] + 4 in our analysis of nonlinear stability. 

(HI) The eigenvalues of A\± are (i) distinct from the shock speed s = 0; (ii) of common 
sign; and (iii) of constant multiplicity with respect to U. 

(H2) det(dF 1 ([/±)) / 0. 

(H3) Local to U(-), stationary solutions of (II. ID . connecting £/±, form a smooth manifold 
{^(•)},(5eWcM ! . 

(H4) The eigenvalues of ^ • £,jdFi{U±) have constant multiplicity with respect to £ G M d , 

Structural assumptions (A1)-(A3) and (H0)-(H2) are satisfied for gas dynamics and MHD; 
see discussions in |MaZ4l IZ3| IGMWZ51 IGMWZ6j . 

Alternative Hypothesis H4'. The constant multiplicity condition in Hypothesis (H4) 
holds for the compressible Navier Stokes equations whenever is hyperbolic. However, the 
condition fails always for the equations of viscous MHD. In the paper, we are able to treat 
symmetric systems like the viscous MHD under the following relaxed hypothesis. 

(H4') The eigenvalues of ^ ■ ^jdF^{U±) are either semisimple and of constant multiplicity 
or totally nonglancing in the sense of [GMWZ6J, Definition 4.3. 

Remark 1.1. There will be easily seen that our results also apply to the case where the 
characteristic roots satisfy a (BS) condition^ (see Definition 4.9, [GMW Z6]). a more general 
situation than the constant multiplicity condition, ensuring that a suitable generalized block 
structure condition is satisfied. See Remark 12.31 for further discussion. 

Remark 1.2. Here we stress that we are able to drop the following structural assumption, 
which is needed for the analyses of |Z21 E31 IZllGMWZl) . 

(H5) The set of branch points of the eigenvalues of (A 1 )" 1 (it A + X^yi ^j^)±i ^ G M, 
£ G is the (possibly intersecting) union of finitely many smooth curves r = T]q(£), on 
which the branching eigenvalue has constant multiplicity s q (by definition > 2). 

1.2. Shock profiles. We recall the following classification of shock profiles. 

Hyperbolic Classification. Let i+ denote the dimension of the stable subspace of 
dF l (U+), i- denote the dimension of the unstable subspace of di ?1 (L r _), and i := i+ + 
Indices i± count the number of incoming characteristics from the right/left of the shock, 
while i counts the total number of incoming characteristics toward the shock. Then, the 
hyperbolic classification of profile U(-), i.e., the classification of the associated hyperbolic 
shock (U-,U+), is 

Lax type if i = n + 1, 

Under compressive if i < n, 
Overcompressive if i > n + 2. 



Thanks to one of the referees for his pointing out this extension. 
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In case all characteristics are incoming on one side, i.e. i + = n or i_ = n, a shock is 
called extreme. 

Viscous Classification. A complete description of the viscous connection requires the 
further compressibility index I, where I is defined as in (H3). In case the connection is 
"maximally" transverse: 

1 Lax or undercompressive case 

i — n overcompressive case 

we call the shock "pure" type, and classify it according to its hyperbolic type. Otherwise, 
we call it "mixed" under /overcompressive type. Throughout this paper, we assume all vis- 
cous profiles are of pure, hyperbolic type. 

For further discussions, see \Z2\ Section 1.2] or [Z3, Section 1.2], and the references 
therein. 

1.3. The uniform Evans stability condition. The linearized equations of (11. lj) about 
U are 

(1-7) U t = LU := J2(B jk U Xk ) Xj - Ytf?V) ai 

j,k j 

with initial data U(0) = U . Here, B^ k := B^ k {U{xi)) and AW := dF^(U( Xl ))U - 
[dB3\U{ Xl ))U]U n {xi). 

A necessary condition for linearized stability is weak spectral stability, defined as nonexis- 
tence of unstable spectra 5RA > of the linearized operator L about the wave. As described 
in |Z2}IZ3]. this is equivalent to nonvanishing for all £ 6 M'' -1 , KA > of the Evans function 

D L (l A), 

(see equation (1A.3D in Appendix [A]) a Wronskian associated with the Fourier-transformed 

eigenvalue ODE. Let £ = (£, A). Introduce polar coordinates ( = p(, with ( = (f, A) G S . 
We also define = S d f]{^e\ > 0}. 

Definition 1.3. We define strong spectral stability as uniform Evans stability: 

(D) Dl(^,p) vanishes to precisely I th order at p = for all £ € Si and has no other zeros 
in Si x where I is the compressibility index defined as in (H3) and (|1.6f) . 

The spectral stability of arbitrary-amplitude shocks can be checked efficiently by numer- 
ical Evans computations as in |HLyZl[ HLyZ2| . 

1.4. The GMWZ result. We recall the recent result of Gues, Metivier, Williams, and 
Zumbrun for low-frequency regimes, and refer the reader to their original papers for the 
detail of statements and the proof. 

Theorem 1.4 ( [GMWZ6] . Theorems 3.7 and 3.9; [GMWZlj . Section 8). Assume (Al)- 
(A3), (H0)-(H3), and (H4'). 

Then, the strong spectral stability condition (D) implies the L 2 uniform stability estimate 
for low-frequency regimes (precisely stated below, (|2.13j) . Section W. 
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Example 1.5 ([GMWZ6J, Section 8). Fast Lax' shocks for viscous MHD equations satisfy 
the structural assumptions of Theorem 11.41 

However, it is also shown that 

Counterexample 1.6 ( [GMWZ6] . Section 8). Slow Lax' shocks for viscous MHD equa- 
tions do not satisfy the structural assumption (H4'), and thus Theorem 11.41 does not apply 
to these cases. 

1.5. Main results. Our main results are as follows. 

Theorem 1.7 (Linearized stability) . Assuming (Al)-(AS), (HO)-(HS), (H4'), and (D), we 
obtain the asymptotic L 1 n i/[( rf-1 )/ 2 ]+ 2 — > LP stability of (|1.7|) for all three types of shocks 
in dimensions d > 3, for any 2 < p < oo, with rates of decay 



Theorem 1.8 (Nonlinear stability). Assuming (Al)-(AS), (HO)-(HS), (H4'), and (D), we 
obtain the asymptotic L 1 H H s — > LP n H s stability for Lax or overcompressive shocks in 
dimension d > 3 and undercompressive shocks in dimensions d > 5, for s > s(d) as defined 
in (HO), and any 2 < p < oo, with rates of decay 



provided that the initial perturbations Uq := Uq — U are sufficiently small in L n H s . 

Remark 1.9. The price of dropping Hypothesis (H5) is that the obtained rate of decay 
is degraded by i 1 / 4 as comparing to those established in \L2\ IZ3|, IZ4] or Theorem 11.101 
below. Therefore the rates are possibly not sharp. In fact, we believe that the sharp rate 
of decay in L? is rather that of a d-dimensional heat kernel and the sharp rate of decay in 
L°° dependent on the characteristic structure of the associated inviscid equations, as in the 
constant-coefficient case [HoZl, HoZ2 . 

Our next main result addresses the stability for the two-dimensional case that is not 
covered by the above theorems. We remark here that as shown in |Z3j . page 321, Hypothesis 
(H5) is automatically satisfied in dimensions d = 1, 2 and in any dimension for rotationally 
invariant problems. Thus, in treating the two-dimensional case, we assume this hypothesis 
without making any further restriction on structure of the systems. Also since it turns out 
that the proof does not depend on the dimensions, we state (and prove) the theorem in a 
general form as follows, recovering previous results of K. Zumbrun (see |Z3j Theorem 5.5]) 
for "uniformly inviscid stable" Lax or over-compressive shocks with same decay rates. 




(1.9) 



\u(t) - u\ LP < c(i + ^-W-VjO+Iht-j 

\U(t) - U\ H s < C(l + t)-^\U \ L i nH s, 
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Theorem 1.10 (Two-dimensional case or cases with (H5)). Assume the same hypotheses as 
in Theorems \1.7\ and \1.8\ with additional assumption (H5). Then Lax or over-compressive 
shocks are asymptotically nonlinearly L 1 n H s — > IP n H s stable in dimensions d > 2, for 
any 2 < p < oo, with rates of decay 

\u(t) -u\ LP < c(i + tr^-^m^s 

\U(t)-U\ H s<C{l + ty^r\U \ L , nH s, 



provided that the initial perturbations Uq := Uq — U are sufficiently small in L PiH s . Similar 
statement can be stated for linearized stability with same decay rates. 

1.6. Discussion and open problems. As observed in \Z3\ IZ4| . the high-frequency esti- 
mate on the solution operator has already been established without the structural assump- 
tions (H4)-(H5), mainly relying on the damping energy estimates. Hence we shall use it here 
as a black box. We would like to draw the reader's attention to our recent work in [NZ2] for 
a great simplification of this original high-frequency argument, requiring higher regularity 
of the forcing / (to credit, the simplification was based on an argument introduced in [KZJ 
for relaxation shocks). 

The difficulty of relaxing Hypothesis (H4) and dropping (H5), extending results in \Z2\ 
\Z3\ IZ4] obtained by pointwise bound approach, is that there and in [GMWZl] the authors 
apply the diagonalization of glancing blocks, where the hypotheses are required, to obtain 
rather sharp bounds on resolvent kernel and resolvent solution. We rather use the L 2 
stability bound more directly, avoiding to get sharp bounds on the adjoint problem where 
the diagonalization of glancing blocks must be applied (see Section 12, [GMWZi] ). and as a 
consequence, avoiding the diagonalization error (denoted by j3 in [GMWZlJ or 72 in |Z3| ) at 
the expense of slightly degraded decay, comparing to those reported in \L2\ |Z3, GM WZI] . 
However, the loss i 1//4 of decay is still sufficient to close our analysis for dimensions d > 3 
in the Lax or overcompressive case and for d > 5 in the undercompressive case. As already 
mentioned at the beginning of the paper, this L 1 — > LP resolvent bound will be the key to 
the improvement. 

Our analysis indeed applies to all applications covered by the GMWZ small viscosity 
theory. Hence, the remaining open problem is to treat cases that are not covered by the 
GMWZ theory, that is, the cases when the structural assumption (H4') of Theorem ll.4l is not 
satisfied or more generally when the generalized block structure fails. Counterexample 11.61 
is showing one of such interesting but untreated cases, violating the structural assumption 
(H4') ; 

It is also worth mentioning that the undercompressive shock analysis was carried out 
in |Z3| only in nonphysical dimensions d > 4, and thus still remains open in dimensions 
for d < 3 for systems with or without assumptions (H4)-(H5). Finally, in our forthcoming 
paper |N2| . we have been able to carry out the analysis for boundary layers in dimensions 
d > 2, extending our recent results in [NZ2j to systems with variable multiplicities. It turns 
out that the analysis for the boundary layer case is quite more delicate than those for the 
case of Lax or overcompressive shocks that we are studying here. 
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2. Linearized estimates 
The linearized equations of (jl.ip about the profile U are 

(2-1) U t = LU := J2(Z jk U Xk ) Xj - J2( AiU U 

j,k j 

with initial data U(0) = Uq. 

Then, we obtain the following proposition. 

Proposition 2.1. Under the hypotheses of Theorems \1.7\ and \l.S\ the solution operator 
S(t) := e Lt of the linearized equations may be decomposed into low frequency and high 
frequency parts (defined precisely below) as S(t) = Si(t) + S2(t) satisfying 



(2.2) I^C^C^I/Ls <C(1 + t)-^( 1 -Vp)4-|-^-(i-«)^| / |^ 

for all 2 < p < oo, d > 3, and (3 = {(3\,(3) with j3i = 0, 1 and a defined as 

,^ „x f for Lax or overcompressive case, 

'II for undercompressive case, 

and 

(2-4) \d2ldlS 2 (t)f\ L2 <Ce-^\f\ HM+m , 

for 7 = (71, 7) with 7! = 0, 1. 

Here, we use the same decomposition of solution operator S(t) as in the article of K. 
Zumbrun [Z3]; see (5.152)-(5.153) in [Z3] or (^32|) below. 



2.1. High— frequency estimate. We observe that our relaxed Hypothesis (H4') and the 
dropped Hypothesis (H5) only play a role in low-frequency regimes. Thus, in course of 
obtaining the high-frequency estimate (|2.4p . we make here the same assumptions as were 
made in [Z3J, and therefore the same estimate remains valid as claimed in (|2.4p under our 
current assumptions. We omit to repeat its proof here, and refer the reader to the article 
|Z3j . (5.16), Proposition 5.7, for the original proof. See also a great simplification in [NZ2J, 
Proposition 3.6 in treating the boundary layer case. 

In the remaining of this section, we shall focus on proving the bounds on low-frequency 
part Si(t) of the linearized solution operator. 

Taking the Fourier transform in x := (x2, ■ ■ ■ , Xd) of linearized equation (12. ip . we obtain 
a family of eigenvalue ODE 

L U 

XU = Lf := (BuU'Y - (AiU)' -i ^ A^U + i ^ B n ^U' 

(2.5) " m m 
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2.2. 1? stability estimate for low frequencies. We briefly recall the procedure (see 
[GMWZ1], page 75-85) of reducing the eigenvalue equations to the block structure equations 
and stating the L 2 estimate for low-frequency regimes by the construction of degenerate 
symmetrizers. 

Let U = (u 1 ,u n ) T a solution of eigenvalue equations, that is, (Lg — X)U = f where 
is defined as in (|2.5p . Following [Z3, Section 2.4], consider the variable W as usual 



,i 



W 



with z := bj^u^. + b^u 1 ^. Then we can write equations of W as a first order system 
(2.6) d xl W = g(x 1 ,X,i)W + F, 

with F := (A^f, 0, f n ) tr , where / = (f,f n ) tr and A* := A\ x - A^ 1 )" 1 ^ 1 ; thus, in 
particular, \F\ < C|/|, for some constant C . It is not necessary for us to carry out in detail 
the form of Q; though, see equation (2.65) of [Z3]. Indeed, we are only interested in the fact 
that bounds in L p of W will give those of U in the same norm. 

We go further as in [GMWZT1 page 75] to write this (n + r) x (n + r) system on R as an 
equivalent 2(n + r) x 2(n + r) "doubled" boundary problem on x\ > 0: 



(2.7) 
where 

(2.8) 



d xl w = g( Xl ,\,£)w + F 

TW = on xi = 

W(xi,\,i) = QY+,W-), 
G(x 1 ,\,i)=( G + _° g ), 

TW = W+ - W- 



with F±(x x ) := F(±xi). 

For small or bounded frequencies (A,^), we use the known MZ conjugation; see, for 
example, [MeZlj or [GMWZl] Lemma 5.1]. That is, given any (A, |) G R d+1 , there is a 
smooth invertible matrix <&(xi, A,£) for x% > and (A,^) in a small neighborhood of (A, £), 
such that (|2.7p is equivalent to 

(2.9) a^y = g+(A,oy + F, f(A,e> = o 

where Q + {\ £) := §{+oo, A, £), W = $Y, F = ^F and FY := T<Z>Y. 
Next, there are smooth matrices V(A, ^) such that 



(2.io) v- l g+v 



H 
P 
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with blocks H(\,^) and 

(2-11) P(\,£)=( P + £ 

satisfying the eigenvalues p of P± in {±5Re/x > c > 0} and 

H(X,i) = H (X,i) + O(p 2 ) 

d 

H (X,i) : = -(^^((iT + 7 )< + £>i4 

i=2 

Define variables Z = (uh,up) T with up := (np + ,up_) T as 

w = $y = f z := r$yz, 

and (F H ,F P ) T = V~ X F. We have 

Let (•, •) denote the standard L 2 product over [0, oo), that is, 

(f,<j) = / /(xi)<7(xi)<ixi, V G L 2 (0,oo), 

JO 

where 3 is the complex conjugate of g. 

Then, recalling that p = |(£,A)| and 7 = KeA, we obtain the maximal stability estimate 
for the low frequency regimes ( [GMWZ61 IGMWZlj ): 

f2 13) {i + P 2 )\uh\ 2 L 2 + \u P+ \\ 2 +p 2 \u P _\ 2 L2 + \u H (0)\ 2 + \u P+ (0)\ 2 + p 2 \u P _(0)\ 2 

<\(SFp + ,u P+ )\ + \(SFp_,up_)\ + \(SF H ,u H )\ 

where S is the degenerate symmetrizer constructed in |GMWZ1] (see equations (8. 2), (8. 2), 
and (6.18)) as follows 

S P 
S H , 



(2.14) S 
and 

(2.15) S P ■■ 



K Id 

-p 2 

for sufficiently large constant K (independent of small parameter p); here, Id is the identity 
matrix and the two subblocks in Sp have the same sizes as those of P± in (|2.1ip . corre- 
spondingly. Here and throughout the paper, by / < g, we mean f < Cg, for some positive 
constant C independent of p. 

There are two possibly subtle points in quoting (|2. 13|) that we would like to point out, 
namely, (i) the estimate (|2.13p was proved in Section 8, [GMWZ1J, under the assumption 
(H4), but not under the relaxed Hypothesis (H4'), and (ii) the estimate was obtained only 
for the Lax shock case. However, in the first matter, the variable multiplicity assumption 
is only involved in the hyperbolic part (the H block in f)2. 12j) ) and the parabolic blocks 
P± remain the same. Thus, the degenerate Kreiss-type symmetrizers techniques (only 
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involved in the parabolic blocks) introduced in |GMWZ1~] can still be applicable here. For 
the hyperbolic part, we now use the recent construction of Kreiss-type symmetrizers in 
[GMWZ6] that applies to the relaxed Hypothesis (H4'), thus yielding the L? estimate for 
this block. In dealing with the second matter, we recall that a crucial step in the analysis 
of [GMWZl] for the Lax shock case was to proving the "right" degeneracy of the boundary 
operator or Lemma 7.1 in [GMWZl] . connecting with the Evans stability condition (D). We 
then observe that with slight modification of the proof, the lemma remains unchanged for 
the under /over-compressive shock case, yielding the same result. For sake of completeness, 
we shall recall the proof of Lemma 7.1, [GMWZl] . with a straightforward extension to other 
cases than the Lax case in Appendix A. 

In other words, with our above observations, we may use the estimate (|2. 13j) as stated 
under our current assumptions in treating all three types of shocks. In addition, thanks 
to that fact that the symmetrizer S is degenerate with order p 2 in the block P_ (see 
(|2.14p . (|2,15p above), we can further estimate f)2. 13|) as 

(j + P 2 )WH\h + \u P+ \\ 2 +p 2 \u P _\ 2 L2 + |n^(0)| 2 + |u P+ (0)| 2 + p 2 \u P _(0)\ 2 

<(\Fp + \,\u P+ \) + p 2 (\Fp_\,\up_\) + (\F H \,\u H \). 

We note that in a final step in [GMWZl, equation (8.11)], the standard Young's inequality 
was used to absorb all terms of (uu,up) hito the left-hand side, leaving the L 2 norm of F 
alone in the right hand side. For our purpose, we shall keep it as stated in (I2.16j) . 

We remark also that as shown in [GMWZlJ, all of coordinate transformation matrices are 
uniformly bounded. Thus a bound on Z = (un,up) T would yield a corresponding bound 
on the solution U. 

2.3. L 1 — > LP estimates. We establish the L 1 — ► LP resolvent bounds for solutions of eigen- 
value equations (L^ — X)U = f in the low frequency regime; specifically, we are interested 
in regime of parameters restricting to the surface 

(2.17) r ? ~:={A : J?eA = -6>i(|£| 2 + |9mA| 2 )}, 

for 0\ > and |(£,A)| sufficiently small. The curve was introduced in |Z2l equation 
(4.26)]. Introducing T f is in fact regarded as a key to the analysis of long-time stability in 
multidimensions. The main point here is that even though A enters into the stable complex 
half-plane ({9?eA < 0}), remains outside of the essential spectrum of limiting linearized 
operators L^ ± ; see \L2>\ Lemma 2.21]. 

In addition, in a related matter, we would like to recall that the Kreiss' symmetrizers 
constructed by O. Gues, G. Metivier, M. Williams, and K. Zumbrun can be attained in a 
full neighborhood of basepoint (£, A) even for IRA = (see, e.g., Theorem 3.7, [GMWZ6J). 
Thus, the L? estimate (|2~16jl is in fact still valid in any region of 

7>-0(M 2 + ieT) 

for sufficiently small. In particular, we shall use (|2.16p for A restricted on the curve 
We obtain the following: 
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Proposition 2.2 (Low-frequency bounds). Under the hypotheses of Theorem \1.8\, for A € 

and p := A)|, 9\ sufficiently small, there holds the resolvent bound 

(2-18) |(L ? - - Ar^/U^O < Cp-^+^-^l/U!^), 

/or all 2 < p < oo, f3 = 0, 1, and a defined as in (|2.3p . 



Proof. Changing variables as above subsection and taking the inner product of each equation 
in (12.12P against uh and up ± , respectively, and integrating the results over [0, x\], for x\ > 0, 
we obtain 



(2.19) 



I 1 Z^' 1 

-|«F/(a;i)| 2 = -\u H (0)\ 2 + Ke / {H(\,£)u H ■ u„ + F H ■ u H )dz, 

II J o 

1 1 

-|u P± (2;i)[ 2 = -|n P± (0)| 2 + »e / (P±(A,0«P± ■ «p± + ^ -u P± )dz. 

Jo 



This together with use of Young's inequality into the last terms involved in i* 1 and the 
facts that \H\ < Cp and \P±\ < C yields 

^ I«h|1<»( x1 ) i$ l^(0)| 2 + p|l*tf|£a + \F H \\i, 

l«p±li'»(*i) < kp±(°)l 2 + I«p±Il2 + IPp±I|i- 

We are now in position of applying the I? stability estimate (|2.16|) . In (|2.20p . multiplying 
both sides of equations of uh by p, of up + by 1, and of up_ by p 2 , adding up results, and 
applying (|2.16p . we obtain 

P 2 (\ U H\ 2 L 2 + \Up\l 2 ) + p\u H \ 2 L ac + \U P+ \ 2 LX + p 2 \u P _\ 2 L oo 

l2JU Z <l*p + [, l«p + |) +P 2 (|PpJ, l«p_[) + (|Pp|, M) + [FhIL + |Pp||i, 

(noting that /? is assumed to be small; in particular, p < 1.) 



Applying again the standard Young's inequality: 

( | F H | , | u H I ) + ( I F P+ 1 , | u P+ 1 ) + p 2 ( | F P _ | , | u P _ | ) 



< 



+ |wp+||°° +p 2 I«p_Il^ 



+ a 



- l \F H \ 2 +\F P+ \l x +p 2 \F P _\ 2 



with e > being sufficiently small, from (|2.21|) . we easily arrive at 

P 2 (\U H \ 2 L2 + \U P \ 2 L2 ) + p\u H \ 2 Loo + |u P+ ||oo +p 2 \u P _\ 2 L oc 

< P^lFnlh + \F P+ \ 2 L1 + p 2 \F P _\ 2 Ll + \F H \ 2 L1 + \F P \ 2 L1 . 
Therefore in term of Z = (uh,up) , simplifying the above yields 
(2-23) P 2 \Z\h (xi) + p 2 \Z\ 2 Lao(xi) < Cp-Vll* 

Now from the change of variables Z = V~ l $~ l W , we have the same estimates for W and 
thus U, because all coordinate transformation matrices are uniformly bounded. Hence, we 
also obtain bounds (I2.23h for U or by the interpolation inequality: 

(2-24) \U\ LP{xi) <Cp- 3 / 2 \f\ L1(xi) . 

This thus proves the proposition in the case of (3 = 0. 
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For (3 = 1, we expect that d xi f plays a role as u pf" forcing. Recall that the eigenvalue 
equations (L^ — X)U = d Xl f read 

LoU 

(B n U Xl ) Xl - (A l U) Xl -i^A^U + i^B' 1 ^ 

(2.25) m m 

+ iJ2( Blk ^U) Xl - J2 Bjk ^kU -XU = d x J. 

Now modifying the nice argument of Kreiss-Kreiss presented in [KK^ lGMWZl]. we write 
U = V + Ui, where V satisfies 

(2.26) (L - X )V = 8 Xl f, Xl el, 

for Ao = p- Noting that A 1 and B n depend on x\ only, we thus can apply here the one- 
dimensional Green kernel bounds investigated by C. Mascia and K. Zumbrun as follows. 
Let G\ a be the Green kernel of Ao — Lq. Observe that our assumptions as projected 

on one-dimensional situations (i.e., £ = 0) are still the same as those in [Z3| . Thus, we 
apply Proposition 4.22 in [Z3| for (I2.26p . noting that Ao = p is sufficiently small. After a 
simplification, we simply obtain 

(2.27) \Gl { Xl ,y x )\ < CV^e"^^-^ 1 ! + e -p\'^-v\ 
and 

(2.28) \d yi Gl { Xl , yi )\ < C[p^e- e ^(pe- ep M + ae~ e ^) + e- p \ x ^{p + ae" 9 ^)}, 

where a is defined as in (j2.3f) . We would like to remark here that Lemma 5.23, [Z3J, 
gives the estimate (I2.28P with a = only for the Lax or overcompressive shocks. For the 
undercompressive shocks, we must have the weaker bound by the term e _6, ' yi ', that is, a = 1 
(for further discussion, see, e.g., equations below (5.106), [Z3]). 

Hence, using (I2.28P and applying the standard Hausdorff- Young's inequality, we obtain 

(2.29) \V\lp( x1 ) + \V X1 \ LP ( X1 ) < I/Ili^i) + ^ P~ a \f\L 1 {x 1 ), 

for all 1 < p < oo and a = or 1 defined as in (|2.3j) . 

Now from U% = U — V and equations of U and V, we observe that U\ satisfies 

(2.30) {L i -X)U 1 = L(V,V X1 ), 
where 

L(V,V X1 ) : =iJ2A j H j V-iJ2B jl £ j V Xl -i^2(B lk CkV) Xl 

m m k& 

+ E Bjh ^v + (A - x )v 

= pO(\V\ + \V Xl \). 



Id 
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Therefore applying the result which we just proved for j3 = to the equations (|2.30[) . we 
obtain 



U X \ LP{X1) < Cp- 3 / 2 \L(V,V xl )\ L i {xi) < Cp-^p \ V \v + \V xl \v 



(2.31) 



3/2 



< Cp 



-3/2+(l-o)| 



\f\V-{x x )- 

Bounds (|2.29p . (|2.3ip on V and U\ clearly give our claimed bounds on U by triangle 
inequality: 

\U\lp < \V\lp + \Ui\ L p. 
We thus obtain the proposition. □ 

Remark 2.3. Under the general structural assumptions, our proof of the L 1 — > L p bounds 
above depends only on the L 2 maximal estimate (12. 161) . As the GMWZ theory covers to a 
more general case than (H4'), namely, the (BS) condition (Definition 4.9, [GMWZ6]), our 
results thus apply to this case as well without any additional work. 

Remark 2.4. In Appendices [B] and O we will prove a slightly-weaker resolvent estimate 
like (f27[8]) in which the Green kernel bounds ([237]) . ([2728]) will not be used. Thus, our 
main results can in fact be derived completely independent of the pointwise Green function 
estimates. 

2.4. Estimates on the solution operator. In this subsection, we complete the proof of 
Proposition 12.11 As mentioned earlier, it suffices to prove the bounds for S\{t), where the 
low frequency solution operator S\(t) is defined as 

1 



(2.32) 



Si(t) :-- 



(2m) c 



|£|<r -/r«n{|A|<r} 



Proof of bounds on S±(t). We first prove (|2.2p for (3 = 0. Let u(xi,£, A) denote the solution 
of (Z/£ — X)u = f, where f(xx,£,) denotes Fourier transform of /, and 



u(x,t) :=5i(t)/ 



(2m) d y^[< r /r«n{|A|<r} 



e^^L^-Xy'fix^OdXdl 



Using Parseval's identity, Fubini's theorem, the triangle inequality, and Proposition 12.2 
we may estimate 



u 



L 2 (xi,x) 



1 



< 



(2tt) m 
1 



xi Jm<r 



e xt u(xi,^,X)dX d£dxi 



r«n{|A|<r} 



MeAt I 



r«n{|A|<r} 



^ xt p^/ 2 dX 



(2TT) 2d 
< C\f\l Hx) 

Specifically, parametrizing by 

X(i,k) = ik-e 1 (k 2 + \i\ 2 ), ke 



xi,i,X)\L 2 (xi)d^ d£ 



r«n{|A|<r} 



2 _ 
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and noting that \d\/dk\ is bounded on n {|A| < r}, we estimate 



15 



e XeM p -3/2 dx 



r«n{|A|<r} 



2 d£ < C 



<C y_ e -2%ICl^|||-l-2e 

noting that f R d-i e~ 9 ^ 2 \x\~ a dx is finite, provided a < d — 1. 
Similarly, parametrizing as above, we estimate 



-0i& 2 t|j,|e-l 



dk 



di 



£ Jr«n{|A|<r} 



< C|/ln ( *) / i . e^ xt p-V 2 d\di 
Jr^n{\\\<r} 



d-l | 1 
< Ct _ — + 4. 



€| 2 t|?|-l/2- e 



k\ e ~ dkd£ 



The xi— derivative bounds follow similarly by using the version of the L 1 — > LP estimates 
for [3\ = 1, noting that in the undercompressive case, both f3\ = and f3\ = 1 have the same 

bounds. The x— derivative bounds are straightforward by the fact that d~f = (i£)^f- □ 



2.5. Proof of linearized stability. Applying estimates (|2.2p and (|2.4p on low- and high- 
frequency operators Si(t) and 52 (t) obtained in Proposition 12.11 we obtain 



(2.33) 



\U(t)\ L 2 < IS^UqIli + \S 2 (t)U \ L 2 

<C(l + t)-^\U \ L i +Ce- r ' t \U \ L 2 
<C(l + t)-^\U \ L inLi 

< C\f\ H s/ K d\ for s > d/2; see, for 



and (together with the Sobolev embedding: \f\ L 
example, [241 Lemma 1.4]) 



(2.34) 



\U(t)\ L oo < \Sx(t)U \ L oo + \S 2 (t)U \L°° 

"^I^oIl 1 + C\S2(t)Uo\ H [(d-i)/2]+2 

^I^oIl 1 + Ce~ vt \U \ H [(d-i)/2\+2 



d-l i i 

< c{i + t)-— 

< c(i+ty— 

< C(l + t)"T- + i|l7o| 



These prove the bounds as stated in the theorem for p = 2 and p = oo. For 2 < p < oo, we 
use the interpolation inequality between I? and L°° . 
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3. Nonlinear stability 

Defining the perturbation variable U := U — U, we obtain the nonlinear perturbation 
equations 

(3.1) U t -LU = J2Q J (U,U x ) Xj , 

3 

where 

Q j (U,U x ) = 0(\U\\U x \ + \U\ 2 ) 
Q J (U,U X ) X] = o(\u\\u x \ + \U\\U XX \ + \U X \ 2 ) 
so long as \U\ remains bounded. 

Proof of Theorem \1.8[ We prove the theorem for the Lax or overcompressive case. The 
undercompressive case follows very similarly. Define 

(3.3) C(t)-= sup (|[/(s)| L2 (l + s)^ + \U( a )\ L ~(l + s)^- l A. 

We shall prove here that for all t > for which a solution exists with ((t) uniformly 
bounded by some fixed, sufficiently small constant, there holds 

(3-4) C(t)<C(\U \ LlnH s+at) 2 ). 
This bound together with continuity of C(i) implies that 

(3.5) C(t) < 2C\U \ L i nHs 

for t > 0, provided that \Uo\ii n us < 1/4C 2 . This would complete the proof of the bounds 
as claimed in the theorem, and thus give the main theorem. 

By standard short-time theory/local well-posedness in H s , and the standard principle 
of continuation, there exists a solution U S H s on the open time- interval for which | C7" | 
remains bounded, and on this interval Q{t) is well-defined and continuous. Now, let [0, T) be 
the maximal interval on which \U\h s remains strictly bounded by some fixed, sufficiently 
small constant 5 > 0. By an auxiliary energy estimate in \Z3\ Proposition 5.9] and the 
Sobolev embeding inequality |C/"|pj/2,oo < C\U\h s (again, see for example, \Z4j, Lemma 1.4]), 
we have 

\U(t)\% s <Ce- et \U \% B + C f e- e ^\U(r)\l 2 dr 

(3.6) Jo 

<c(\u \ 2 Ha +at) 2 )a+tr {d - 2)/2 . 

and so the solution continues so long as ( remains small, with bound (|3.5p . yielding existence 
and the claimed bounds. 

Thus, it remains to prove the claim (|3,4p . By Duhamel formula 

(3.7) U(x,t) =S(t)U + f^S{t-s)Y,dx 3 Q J {U,U x )ds, 
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where U(x,0) = Uq(x), we obtain 

(3-8) \U(t)\ L 2 <\S(t)U \ L 2 + f |5i(t - s)d Xj Qi(s)\ L 2ds + f \S 2 (t - s)d Xj Q? {a)\vda 

Jo Jo 



where \S(t)Uo\i2 < C(l + t) 4 |[/o|LinL 2 as i n the proof of linearized stability, 
ft ft 

/ \S l (t-s)d Xj Q^s)\ L 2ds<C {l + t-8)-*?-*\Q?(a)\ L xda 
Jo Jo 

<C [ {l+t- s)-^r-^\U\ 2 m ds 
Jo 



<C(\U \ 2 H s+((t) 2 ) I (l + t-s)-^-5(l + a) 
io 

<C(l + t)-^(\Uo\ 2 H s+((t) 2 ), 



and 

/ \S 2 (t - s)d x .Q3(s)\ L2 ds < f e- 6 ^\d x .Q\s)\ L2 ds 
Jo Jo 

< C f e- e ^\U\ 2 H sds 

Jo 

< C(\U \ 2 H s +((t) 2 ) fe-^d + sy^ds 

Jo 

<c(i + t)-^(\u \ 2 H s+c(t) 2 )- 

Thus, dividing by (1 + t) ~ , we obtain 

(3.9) \U(t)\ L2 (l + t)^ <C(\U \ L i nHs +C(t) 2 ). 

Similarly, we estimate the L°° norm of U. By Duhamel's formula (13. 7p . we obtain 



(3.10) 



\U{t)\ L oo <\S{t)U \ L oo + j \Siit - s)d Xj Q j {s)\ L °°ds 
+ J \S 2 {t-s)d X] Q j {s)\ L oods 
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where \S(t)U \L°° < C(l + 1) \U \ L i nH [(d-i)/2]+2, 
J \Si(t - s)d X] Q j (s)\ L ocds 

Jo 

<C f\l + t-s)-^ +l ^\U\ 2 Hl 
Jo 

< C(\U \h + C(t) 2 ) f\l + t - s)-^-Hl + s)-*? 
Jo 

<c(i+tr^M(i£/oi^ + c(t) 2 ) 

and (by the Moser inequality; see, for example, inequality (1.22), |Z4|), 

J \S 2 (t - s)d X] Q j (s)\ L oods 
t 

\S 2 (t - s)d Xj Q j (s)\ H[{d -i)/2]+2ds 

< f e- e ^\d x Q\s)\ H[{d -, )m+2 ds 
Jo 

<C [ e- fl (*-*)|l7Uco|l7| H[( a-i )/2]+ 4ds 
Jo 

< C(\U \ 2 H s + C(t) 2 ) f e- e ^ s \l + s)-^ + Hl + sy^ds 

Jo 

<C(l + t)-^ + *(\Uo\ 2 HS + ((t) 2 ). 
Therefore we have obtained 
(3.11) \U(t)\ L ~(l + t) d -^-i <C(\U \ L inHs+((tf) 

and thus completed the proof of claim (|3.4|) . and the theorem. □ 

4. Two-dimensional case or cases with (H5) 

In this section, we give a proof of Theorem ll.101 Again, notice that the only assumption 
we make here that differs from those in [Z3] is the relaxed Hypothesis (H4'), treating the 
additional case of totally nonglancing characteristic roots, which is only involved in low- 
frequency estimates. That is to say, we only need to establish the L l — > L p bounds in 
low-frequency regimes for this new case. We give the proof of these bounds by modifying 
the proof in [GMWZ1, Section 12] and thus will not cite the estimate (|2.16p in this section; 
in fact, the proof is completely independent of previous sections. In addition, our proof is 
somewhat more direct and simpler than those in [GMWZ1, Section 12] by not bypassing to 
the dual problem. 
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Proposition 4.1 (Low-frequency bounds; |Z3j . Corollary 5.11). Under the hypotheses of 
Theorem \1.1(A for A £ (see (|2.17p ) and p := A)|, 9\ sufficiently small, there holds 
the resolvent bound 

(4-1) |(% - \y l d^f\ LHxi) < CW^I/Il^O' 

/or all2 < p < oo, (3 = 0, 1, and 72 is i/ie diagonalization error (see [Z3| . (5.40)) defined as 

(4.2) 72 := l + £ [p-^mA- r?±(C)| + p] 1/Sj_ \ 

j'.± 

We again perform the standard procedure (see Section I2.2[) of writing the linearized 
equations in form of the first order eigenvalue equations f)2 . 12|) : 

Locally, in a neighborhood of a base point Xq := 0) with £ = (r, 7, £) and A = 7 + it, 
we further use the Assumption (H4') to write H in block-diagonal structure (see [GMWZ1, 
Proposition 6.1]) with appearance of a new mode, totally nonglancing, and decompose the 
resolvent solution U into 

(4.4) U = Up + U He + U Hh + U Hg + U Ht , 

corresponding to parabolic, elliptic, hyperbolic, glancing, or totally nonglancing blocks. We 
further write 

U t = U i+ + Ui- 

for i = P, H e , Hh, H g , Ht, where Ui± are defined as the projections of Ui onto the growing 
(resp. decaying) eigenspaces of Q + in (|2.10p with respect to the corresponding blocks. 

These first four blocks have been treated in |Z3} Corollary 5.11] or [GM WZlj Corollary 
12.2] for which the totally nonglancing modes are absent. For sake of completeness, we 
treat these modes again here in a slightly different analysis, modifying those of [GMWZ1, 
Section 12]. In fact, since each mode interacts with the other via the Evans condition (D) 
or, more precisely, the boundary estimate (|A.28j) . we cannot obtain (|4.ip for each mode 
separately. 

We shall use the following simple lemma. 

Lemma 4.2. Let U be a solution of d z U = QU + F with U(+oo) = 0. Assume that there 
is a positive [resp., negative] symmetric matrix S such that 

(4.5) "51SQ := ^(SQ + Q*S*) > Old 

for some 9 > 0, and S > Id [resp., —S > Id]. Then there holds 

\U\ho+0\U\b<\F& 
[resp., \U\ 2 Lao +e\U\ 2 L2 < \U(0)\ 2 + \F\ 2 L i J. 
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Proof. Taking the real part of the inner product of the equation of U against SU and 
integrating the result over [xi,oo] for the first case [resp., [0, x\] for the second case], we 
easily obtain the lemma. □ 

Proof of Proposition \4~J\ As in [GMWZl, Section 12.2], the first step is to put blocks into 
a diagonal form; indeed, parabolic blocks are already diagonal as in (j2.11[) : hyperbolic 
blocks are lxl blocks with real part vanishing at the base point Xq, but with real part 
> (resp. < 0) when p > in polar coordinates (thus, vanishing at order p 2 in original 
coordinates); elliptic blocks are those Q k with "$.eQ k positive or negative definite at the base 
point (thus, vanishing at order p); and finally glancing blocks are of size larger than lxl 
whose components are purely imaginary at the base point. We recall the following lemma 
in [GMWZ1], diagonalizing these glancing blocks. 

Lemma 4.3 (Lemma 12.1, [GM WZl] ). Diagonalize the glancing blocks Q k by the transfor- 
mation Tjj , where Th may be chosen so that 



(4.7) \T Hg \<C, |T^|<C 72 , |T-; |ir _|<(7 7 



(4-10) 



n 

where 72 is defined as in (I4.2p and 71 is defined as 

. . r , + .~. i (i-[(u k +i)/2])/u k 

(4.8) 7 i := max p \ismX - rj, (£)| + p , 

k L 

and T^j ^ denotes the restriction of to subspace H g ^. 
In addition, after a further transformation if necessary, 

(4.9) Q' k := T^Q k T Hg = diag(a fej i, • • • , a k ,i, a k ,i+i, ■ ■ ■ ,a k , Uk ) 
with 

-7f 2 Ke a k j > Cp 2 , j = 1, I, 
7f 2 Ke a kJ > Cp 2 , j = l + l,...,u k . 

Remark 4.4. 71,72 are identical to a, (5 in [GMWZl], respectively, and (|4.10p was calcu- 
lated in |GMWZll equation (12.40)]. 

We now can work in diagonalized coordinates: 

U' := T„ X U 

where Tjj are obtained as in Lemma 14.31 for glancing blocks and identity matrices for 
the other blocks. In these coordinates, since blocks are diagonal and growing/decaying 
subspaces (at least for the first four modes) are separated, we apply Lemma 14.21 for each 
block with S = ±Id, yielding 

\uL\l^ + e l \ul_\ 2 L2 <\uUo)\ 2 + \FU 2 L1 , 

where 6i = 1, p, p 2 , min^ |Ke a k j | for i = P, H e , H^, H g . 

For the totally nonglancing blocks Q\ , as constructed in [GMWZ6], Lemma 5.3, there 
exist symmetrizers S k that are definite positive [resp., negative] when the mode is totally 
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incoming [resp., outgoing]. Denote U' Ht+ [resp., U' Ht _] associated with totally incoming 
[resp. outgoing] modes. Then by applying Lemma 14.21 with 9 = p 2 , we obtain 

\TT> I 2 _l_ r?\TT> I 2 < I p/ |2 
( . , \ U H t+ \L°°+P \ U H t+ \L* ~ \*H t+ \L^ 

{ ' ) \ukJU+p 2 \ukJh<\ukAo)\ 2 + \FkJl- 

To finish the proof, we only need to deal with the boundary terms |C/|_(0)| 2 and \Ujj t _ (0)| 2 
in (|4.1ip . (|4.12p . We could use a more detailed version of the L 2 stability estimate (|2.16p . 
corresponding to each diagonal blocks (see [GMWZ6J), yielding bounds on these boundary 
terms. However, let us now follow the boundary treatment presented in [GMWZ1, Section 
12.3] instead, being rather independent of (|2.16p . 

The diagonalized boundary condition is V := TTh 9 - By computing, we observe that 

\v'u' Hg _\ = \vu Hg _\ > c- l \u Hg _\ > C lT -i Hs ] 1 > c- l ii l \u' Hg _\. 

Thus, together with (1X28]) 

(4.13) \T'UL\ = \TU-\ > C- 1 [\U' He _ | + \U' Hh _ | + \U' Ht _ | + ^U'^ I + pW'p_ 
Meanwhile, we have at x\ = 

(4.14) \T'Ui\ < \T'U'\ + \T'U' + \ < \U' + (0)\ < \U' + \l*,. 

Now, multiplying the first equations in ([4. lip . (|4.12p by a sufficiently large constant k 
and the second equations by l,7f 2 , or p 2 , corresponding to each block with its boundary 
degeneracy of order in (|4.13p . and adding up the results, we easily obtain 

, . x \ u 'p+\l°° + \ u 'p+\l 2 + \ u 'h+\l™ + P 2 \ u h+\l 2 + < If'I 2 

{ ' J2ittI |2 i „2\tt) |2 _i_ n ,-2\TT> |2 , „2\ttI |2 ~l li 1 ' 

P \ u p-\l°°+P I^p-Il2+7i \Uh-\l™+P \ u h~\l2 
(noting that p < 1, p < j^ 1 < 1, and 7-f 2 minj |3f?e ak,j \ ^ P 2 by (|4.1U|) ). This yields 
(4-16) p 2 \Ur L ~+p 2 \U'\l*<\F% 
or equivalently, 

(4.17) \U'\lp < p'^F'Ili, Vp>2. 

Thus, by recalling that U = T Hg U' and F' = T^F, KT7\\ and (gZI) immediately yield the 
proposition for (3 = 0. For = 1, we can follow the Kreiss-Kreiss trick as presented in the 
proof of Proposition 12.2^ thus completing the proof of Proposition 14.11 □ 

Proof of Theorem 1 1.1 01 Proposition 14. II is the Corollary 5.1 in |Z3j with an extension to the 
totally nonglancing cases. Thus, we can now follow word by word the proof in [Z3], yielding 
the theorem. □ 

Remark 4.5. We have seen in the above argument that the existence of positive/negative 
Kreiss' symmetrizers with an appropriate constant 6 (in Lemma l4.2p is sufficient to obtain 
the result. Though, proving the existence of such symmetrizers is a highly-technical task in 
general for variable multiplicity blocks. See [GMWZ5, GMWZ6 . 
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Appendix A. Evans function for the doubled boundary problem 



For sake of completeness, we recall here the proof of Lemma 7.1, [GMWZ1] and its 



(A.l) 



straightforward extension to the case of over- and under-compressive shocks. 
Consider the 2N x 2N doubled boundary problem (|2.7p (with N := n + r) 

U x -G(x,()U = F, 
TU = on x = 0, 

where U = (U+, U-) and TU = U+-U-, with U + = {U 1 ,..., U N ), U- = (U N+1 , U 2N )- 
Let £-((, p), for 7 > and p > 0, be the space of boundary values at x = of decaying 
solutions to the homogeneous problem 

U x -G(x,QU = 0. 

As shown, for example, in |GMWZ6| Theorem 3.7], the space £-((,, p) has a continuous 
extension to a small neighborhood of 7 = 0, p > 0. Then the Evans function for (|A.1[) is 
defined as the 2N x 2N determinant: 

(A.2) D(C,p) = det(kerr,f_)| a=0 . 

Meanwhile, the Evans function Dl for the problem f|2.6j) is defined as 

(A.3) D L ((,p) = det«, . . . MkMk+i, ■ ■ ■ Mn)\x=o 

which is analytic for 3fteA > and can be continuously extended to a small neighborhood 
of SfteA = (see, e.g., Lemma 5.24, |Z3]). Now let <j>j, j = be the derivative of 

the profile U s with respect to Sj, where I is the dimension of the smooth manifold {U s (-)} 
defined as in (H3). Thanks to the Evans condition (D), without loss of generality, we can 
assume that 

(A.4) HjW.O) =Z4_ i+ i(*,C,0) = 0), 

for j = 1,... J. 



Let ej € C be the unit vectors 



N 

(^■(0),0) 



j = 1,... ,1, 



J 10,(0)1 

and extend to an orthonormal basis e%, . . . , e^y of C N . Then the Evans function (|A.2p for 
the doubled boundary value problem can be explicitly defined as 

/ . r \ , . x /ei ... e N UP ... U? ... 

We also set 

(A.6) S-A(,P) = span ( (^fh) (ft 



IlL 1 ■■ • ) l 7/i [ i(o,c,p)" 



For e > fixed, denote by St^^Cp) any complementary subspace in £_(£,p) varying 
continuously with (£, p) such that 

(A.7) e-(S,p)=S-.j(t,p)<B £i,<,, e (C, P) 

with uniformly bounded projections for < p < e. 
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Then, we recall the following proposition that was proved for the Lax shock case in 
[GMWZlj . Proposition 7.1. 

Proposition A.l. (1) Let Dl(C,p) and 1D)(C,, p) be the Evans functions defined as above. 
Then 

(A.8) D L ((,p) = (-l) N m,p). 

(2) Under the Evans assumption (D), we have the following. 

(a) For any choice of < 5 < R there is a constant C$, R such that when 5 < p < R, 

(A.9) \Tu\>C StR \u\ forue£-(C,p). 

(b) There exist positive constants Ci,C2,S such that 

(A.10) Cxp\u\ < \Tu\ < C 2 p\u\ foru G £_^(C,p) 

for < p < 5. 

(c) There exists C > such that 

(A.ll) |ru|>C|u| for u€ St^iC, p) 

for < p < e. 

(d) For any choice of R > there is a constant Cr such that for < p < R, 
(A.12) \Tu\>C RP \u\ forue£-((,p). 

Proof. We follow word by word the proof for the Lax shock case in [GMWZl l] . Proposition 
7.1. First, by performing the row matrix operation, (1) is clear. (2a) follows by continuity 
and compactness, and the fact that Tu is nonvanishing for nonzero u G £-(C, p) when p > 
by Evans function assumption (D). 

For the proof of (2b), let us denote the matrix in (|A,5|) by M and perform column oper- 

/ \ 

ations to replace the last I columns of M. by I L 3 I , and call the resulting matrix 

\ U N-j+lJ 

Now thanks to the normalization (|A.4[) and the fact that fast modes depend analytically on 
p, we have for j = 1, . . . , I 

Thus, the definition of ej, linearity of the determinant in the last I columns, and the Evans 
condition (D) show that cy — c 2 j are nonzero for all j. This together with the definition of 
TU 

P = U* - U L 
yields (|A.10j) at once. 

(2c) Let «i(CjP) 5 ■ • • ,vn{(,p) be the last 2n columns of the matrix A4i defined above. 
These vectors form a basis for f_(^,p). Take an arbitrary vector w G E c _ ^, e (C>p)- Then 

N 

(A- 14) w = J2 c jAC,P)vj{C,p), 

3=1 
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where Cj te ((,p) depend continuously on ((,p). 

Set c e = (ci j£ , • • • , CN-i : e) and c t = (cjv-z+i,e, ■ • • , cjv,e)- The condition that the projec- 
tions in (|A.7p are uniformly bounded implies that there is an eo > such that 

(A.15) \c' € {(,p)\>e \c^(,p)l 
for < p < e. 

In view of (D), we just need to show that Tw is nonvanishing at p = for w as in (|A.14j) 
and (1A.15|1 with |(c e ,c e )| = 1, since (1A,11|) then follows by continuity and compactness. 
Suppose Tw = at (£, 0) for some such w. Because of (jA.15j) some Cj >e with j < N — I, say, 
ci j£ satisfies 

(A.16) \ci, e \>co 

for p near 0, and for some cq > 0. Since Tw = at p = and w(C, p) is continuous, we have 
(A.17) w(C, P )= ( a< $)+0(p). 



AO, 

Write Vj = (vj + ,Vj-), use column operations to replace v\ in M\ by w, and call the resulting 
matrix A4 2 - Then A4 2 = 

e x ... e N a{Q + 0{p) v 2+ ... v N ^ l>+ (&(0), 0) + 0{p) ... (0,(0), 0) + 0{p)\ 
ex ... e N a(() + (D(p) v 2 ^ ... v N ^_ (&(0), 0) + Ofj>) ... (<t>i(0), 0) + O(p)] 

(TAT71) implies that | det M 2 (C, p)\ > C\det Mx((, p)\ for some C > uniformly near (CO). 
But 

det M 2 ((, P) = 0{p) l O{p) as p -» 0. 
This contradicts the assumed vanishing of det = det A4x to exactly Z*' 1 order at p = 0. 

(2d) For any fixed (C>/o) ; let u* = ( n+ ^'^ ] g £_(£, p) be an element where the 

\u-{C,P)J 

minimum 

min \Tu\ 

|«|=l,tee£_ (C,p) 

is attained. Write u* = Y^j=\ c j,e(C> P) v j{C, p) an d define c* e , c* e in the same way as above. 
Then, again, the uniform boundedness of the projections in (1A.7D implies that there is an 
eo > such that either 

(A.18) \c*' e (C,p)\>e \c:"(C,p)\ 
or 

(A.19) \c*"e(C,p)\ >eoK'((,p)\ 

for < p < e. Correspondingly, these imply that, without loss of generality, there holds 
either 

(A.20) \c*xJ>co 
or 

(A.21) \c* N J>co 
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for p near 0, and for some Cq > 0. 

In the case that (|A.20|) holds, as above, we perform column operations to replace v% in 
M\ by it*, and call the result M 5 . Then M 5 = 

'e x ... e N u + v 2+ ... v N ^+ (&(0), 0) + O(p) ... (0,(0), 0) + 0(p)\ 
ei ... e N u_ ... v N -i- (^i(0), 0) + O(p) ... (<&(0), 0) + 0(j>)J ' 

Next perform column operations to replace the column u* = by 

it+ — U-\ _ fYu* 

o v o 

Thus, by direct calculations and (|A.20[) . 

|detM 3 | = \Tu*\0(p) 1 > C| detail = C\detM\ > Cp l , 

which gives \Tu*\ > C. 

Similarly, in the case that (|A.20|) holds, replacing vn in A4\ by u*, denoting the resulting 
matrix by M±, and performing column operations as above, we then obtain 

|detM 4 | = MC) - u^l P )\o{p) 1 - 1 = [r^io^)'- 1 . 

This together with | det A4 4 | > C\detMi\ = 0(p l ) by (|AT2~H yields |IV| > Cp. 

Thus, altogether we obtain \Tu\ > Cp\u\, for u € £-((,p) with \u\ = 1, uniformly in p 
near 0. Together with (2a), this proves (2d). □ 

Now let T be the MZ conjugation such that (|A.ip leads to the following constant- 
coefficient system 

U x -G(ocX)U = F, 
TiU = on x = 0, 

where T\ = TT and G has the block form as in (|2.10p . (j2.11j) : 

/p+(C) o o 

(A.23) G(oo,C)= P_(C) 

V H(C,p), 

Thus, we can decompose U € C 2N as follows 
(A.24) U = U P+ + U P _ + U H+ + U H _ , 

and set 

U- = Up_+U h _ 6f_(C,p). 
Define the Z-dimensional subspace Pp 1 _ of Pp_ by 

£-,</> = TEp 1 _ j 

where is defined as in (|A.6|) . and for e > fixed, chose a smoothly varying complemen- 
tary subspace Ep i _ such that 

fA25) E P _ =E Pl _(C,p)®Ep 2 _ it (t,p), 

U P _ = U Pl _ + Up 2 _ e 



(A.22) 
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with uniformly bounded projections for < p < e. Take 
(A.26) £ c _^ e = T{E P2 _ e ®E H _). 

£- Ae i s then a choice that works in (1A.7|) . 

Then the following is an immediate consequence of Proposition IA.11 

Corollary A. 2. There exist positive constants C\, C2 and 5q such that for < p < 80 

(a) C lP \U Pl _ I < |IW A _ I < C 2 p\U Pl _\, 
(A.27) (6) \T!(Uh. +Up 2 _J\ > C 1 (\U H _\ + \Up 2 _J), 

(c) > Cip\U-\, 

where T\ is defined as in (IA,22j) , These estimates hold uniformly near the basepoint Xq = 

(C,o). 

Thus, we obtain the following lemma which is essential for the construction of degenerate 
symmetrizers. 

Lemma A. 3 (Lemma 7.1, [GMWZlJ). There exists a constant 5 > such that for p 
sufficiently small we have 

(A.28) \TU-\>S(\U H _\+p\Up_\) 
uniformly in a neighborhood of the base point Xq = (£, 0). 
Proof. In view of (|A.27p (a), (b), we have 

|IW_| = \TiU H _ + IW Pl _ + TiUp 2 _J 

>C(\U H _\ + \U P2 _J)-Cp\U Pl _\. 

Adding a sufficiently small multiple of this inequality to the inequality (|A.27|) (c) 

liW-i > c P \U-\ = c P (\u H _\ + \u Pl _\ + pp 2 _j), 

we obtain for p small 

|IW_| >6(\U-\+p\U Pl _\ + \Up 2 _J), 
which implies (jA?28l). □ 

Appendix B. Auxiliary problem 
In this section we consider the n x n system on the whole real line R 
(B.l) L V := (B U V X ) X - (A l V) x = f x 

where A^^B 11 are same as in (|2.25p . Let us recall B 11 = K11 ^11^ • We shall derive an 

estimate slightly similar to (|2.29[) by Kreiss-type symmetrizers techniques in the case of 
Lax and overcompressive shocks. This will be done by modifying the proof in [GMWZ1, 
Section 10.2]; though, our purpose is slightly different and we have to treat the degeneracy 
of the viscosity matrix B 11 as comparing to the identity matrix in [GMWZlJ. Specifically, 
we prove the following: 
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Lemma B.l. Let V = (Vi, V2) £ C n-r x C r be a solution of (jB.ip . We prove that there 
exists a constant C > such that 

\v\ L » < c(\f\ L , + 
1 ' j I^Ilp < cd/kx + |/U« + l/.lw) 

/or any 1 < p < 00. 

Proof. We first integrate the equation (|B.1|) . yielding 
(B.3) S 11 ^ - A X V = f. 

Consider the double 2n x 2n boundary problem onx>0 equivalent to (!B.3j) 

BW X - AW = F 
^ BA ^ TW = on {x = 0} 



where, defining <j)±{x) = 4>(±x) for any function <f> defined on 

W[x) ■■ 



A(x) 



V + (x) 
V-(x) 

A\{x) 
-Al(x] 



™ = {-}%, 

TW = V+-V- 

In what follows, we shall keep track of variables W as (W± , W% , W{ ,W%) € C™ _r xCx 
C n ~ r x C r in the obvious way corresponding to matrix blocks as above. Notice also that B 
is degenerate in W^-blocks. 

Let £-(0) be the space of boundary values of decaying solutions of (|B.4p when F = 0. 
Then, we have 

dim£_(0) = i 

where i is defined at the beginning of Section 11.21 On the other hand, ker T has dimension 
n. Thus, Assumption (D) then implies that kerT and £_(0) have an I = i — n dimensional 
intersection spanned by 

(0i(o),0i(o)),-.- ,(<Mo),<Mo)) 

where functions (f>i(0) are defined as in the paragraph just below (|A,3l) . 
We define an augmented boundary condition T with property that 

(B.6) C 2n = kerfe£_(0). 

Since <j)j form a basis of the tangential space of the smooth manifold {U S { )} defined as 
in (H3), without loss of generality, we assume that the j th component of cj)j is not zero. 
Thus, let us define 

(B.7) fW = {W 1 ,--- ,Wi,V+-V-) 
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where W±, • • • , Wj are the first I components of W 6 C 2n . 
Now we consider the system 

BW X - AW = F 

(B.8) 

TW = on {x = 0}. 

Since any solution of (|B.8|) is also a solution of (|B.4p . we only need to give an estimate 
for solutions of (|B.8p . By using the MZ conjugation [MeZl], there is a uniformly bounded 
transformation C such that by setting W = CZ, (|B,8|) gives 



(B.9) 



B(oo)Z x - A{oo)Z = F 
TZ = on {x = 0} 

where T = TC. Now let us define new variable Y as 



f (I, 



Y := QZ, 



with 



Q := 



n—r 





°2 



V 







'-n-l 





-62 



11 



Then Y solves 



(B.10) 



y = 
y = 
fy 
1-1 



viMQ-iy = f 

(^l + 1 ^2 + ' ^1 ' ^2 ) 

(o,y+,o,y 2 -) 

= on {x = 0} 



where T = TCQ~ L and F = FQ 

Now by view of f|B.5[) . (|l,3p . and (H2), eigenvalues /J,j of each block of „4(oo)Q _1 are 
distinct and nonzero. Thus, by performing a further transformation if necessary, we could 
assume that A(oo)Q^ 1 is diagonal. In these diagonalized coordinates, the system (|B.10p 
consists of In "uncoupled" equations: 



(Y2±)a 



Vj±Yi± = Fi± 
Vj± Y 2± = F 2 ± 



where note that Y%± are the projections of y on the growing (resp. decaying) eigenspaces 
of A associated to eigenvalues Uj±. In particular, ±/Xj± > 0. 
From equations for Y±, it is clear that 



(B.ll) 

Meanwhile, Y 2 ± satisfies 
Y2+(x) 
Y 2 -(x) 



\Y1\1* < \h 



LP 



Vp > 1. 



e ^ + (*-y)P 2+ (y) dy, 



eW-*y 2 _( )+ f* eV-( x -riF 2 -(y) dy. 
Jo 
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Thus, this yields 

(B.12) l^<|^, 

\Y 2 -\lp < \F 2 -\lp + 1*2- (0)|, 

Now since l2_(0) G £-(0), by view of (|B.6[) as bounded projections and the fact that all 
our transformations and their inverses are bounded, we must have 

l*2-(o)| < |fy_(o)| < |fy(o)| + |y+(o)| < |y+| L <». 

Altogether, we obtain 
(B.13) \Y\ LP < \F\ L i + |F|loc, Vp>l 

which proves the first bound in (|B.2p . Estimates for derivatives are then immediate by 
differentiating the equations of Y\ in (|B.10p and by solving equations (|B.10p for Y 2x in 
terms of Y and F. Thus, we have proved the lemma as claimed. □ 



Appendix C. Independence of the pointwise Green bounds 

In this section we comment on independence of the pointwise Green function estimates. 
The high-frequency estimate (|2.4p can be derived entirely from auxiliary nonlinear energy 
estimates as done in |Z4j ; see also Proposition 3.6, |NZ2j . for a great simplification. Whereas, 
the independency for the low-frequency estimate (|2.2p can be seen by first proving the 
following slightly-weaker version of Proposition 12.21 independent of the pointwise bounds 
(|2.27p . (|2,28p . A similar version can be done for Proposition 14.11 

Proposition C.l (Low- frequency bounds). Under the hypotheses of Theorem 1 1 . 8[ for A € 
and p := |(£, A)|, 9\ sufficiently small, there holds the resolvent bound 

(C.l) |(L f - A)-X/| LP(a;i) < Cp-W-^l/U^) + \d Xl f\ LHxi) ), 

for all 2 < p < oo, (3 = 0, \, and a defined as in (|2.3p . 

Proof. Certainly by Proposition 12.21 we only need to prove the bound in the case (3=1. In 
the case of undercompressive shocks, the bound is clear by applying (|C.ip with (3 = 0, a = 1 
and / replaced by d Xl f: 

(C2) - \r l d x J\ LP{xi) < Cp- 3 / 2 \d Xl f\ LHxi) . 

Now for the case of Lax or over compressive shocks, we use the Kreiss-Kreiss trick as 
in the proof of Proposition 12.21 that is, write U = V + XJ\ where V solves the auxiliary 
problem (jB.ip . From the estimate (|B.2p of V and the inequality: \f\L°° < (with 
/(+oo) = 0), we have 

(C.3) \V\ LP + \V x \ L i < C(\f\ L i {xi) + \d x J\ L i (xi) ), Vp > 1. 



Thus, replacing (|2.29p by this inequality and following the proof of Proposition 12.21 we 
obtain the bound ([(XTj) . □ 



30 



TOAN NGUYEN 



Proof of Theorem \1.8\ provided (|C.ip . The resolvent estimate (|C.ip is only weaker than 
(|2.18j) by a stronger norm on /. We thus can certainly follow the proof in Section 12.41 
yielding a low-frequency estimate like (|2,2p . but again weaker by a stronger norm on /, 
namely, |/|x,i( x ) + \d xl f\i J 1 (x)- With this slightly weaker estimate for S±, we can follow word 
by word the proof of the theorem in Section [21 noting that the higher derivatives of / (then 
of U) can then be estimated by the energy estimate (I3.6p . Thus, we obtain the theorem 
without requiring any further regularity on the structures of the system. □ 
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